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1. INTRODUCTION 



Transverse Single Spin Asymmetry (SSA) arises in the scattering of a transversely polarized 
proton off an unpolarised hadron or nucleon if the scattering cross section depends on the direction 
of polarization. The Single Spin Asymmetry for inclusive process + S — )• C + X is defined as 

Ja^ - da^ 

A/v = , . . (1.1) 
da^ + da^ 

where do^^^^ denotes the cross section for scattering of a transversely polarized hadron A off 
an unpolarized hadron B, with A upwards (downwards) transversely polarized w.r.t. the produc- 
tion plane. SSAs significantly different from zero have been observed over last 35 years starting 
with pion production in scattering of polarized protons off unpolarised proton target Large 
SSAs have been measured in pion production at Fermilab ^ as well as at BNL-RHIC in pp^ 
collisions [||]. SSAs have also been observed by the HERMES [|] and COMR^SS [|] collabora- 
tions, in polarized semi-inclusive deep inelastic scattering (SIDIS). The magnitude of the observed 
asymmetries has been found to be larger than what is predicted by perturbative quantum chiomo- 
dynamics (pQCD) [|]. 

Theoretically there are two major approaches to explain the SSA's. One is the twist three 
approach and other is the transverse momentum dependent (TMD) approach which we have used 
in the present work. The TMD approach is based on a pQCD factorization scheme in which spin 
and intrinsic transverse momentum effects are included in parton distribution functions (pdf 's) and 
fragmentation functions (ff 's). One of the difficulties in getting information about the spin and 
transverse momentum dependent pdf 's and ff 's is that very often two or more of these functions 
contribute to the same physical observable making it difficult to estimate each single one separately. 

The study of spin asymmetries requires extension of TMD factorization scheme to polarized 
case. Sivers in early 90's proposed that there exists a correlation between the azimuthal distribution 
of an unpolarized parton and spin of its parent hadron [^. Number density of partons inside proton 
with transverse polarization S, three momentum p and intrinsic transverse momentum of partons 
is expressed in terms of Sivers function A^/^/^t {x, k± ) 

f„/pt{x,k^,S)=fa/p{x,k^) + ^A^f^,/pt{x,k^)S-(pxk^) (1.2) 

S • (p X kj^) gives the correlation between the spin of the proton and intrinsic transverse momentum 
of the unpolarised quarks and gluons. There have been studies on the quark Sivers function in 
SIDIS and the gluon Sivers function in the process p^p — )• DX [||, ||]. In this work, we propose 
charmonium electroduction as another probe of the gluon Sivers function. 



2. FORMALISM FOR ASYMMETRY IN 7/1// PRODUCTION 

We have estimated SSA in photoproduction of chaimonium in the process e + p'^ ^ e+J /y + 
X. At leading order (LO), there is contribution only from a subprocess yg — )• cc. In addition, since 
we are using color evaporation model (CEM) for charmonium production, only one pdf is involved. 
Thus the process under consideration can be used as a clean probe of the gluon Sivers function. 
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Also since the charmonium production mechanism can have implications for this SSA, its study 
can help throw some light on the production mechanism of charmonium as well. 

Charmonium production process can be understood in terms of two distinct steps- production 
of a cc pair (a short distance process) and a subsequent binding of this pair in to charmonium (a 
long distance process). Various methods to describe this non-perturbative evolution of the cc pair 
into charmonium lead to different models of charmonium production. As a first step in our inves- 
tigations of SSA in charmonium production, we have used the Color Evaporation model (CEM) of 
charmonium production. According to CEM, the cross section for charmonium production is pro- 
portional to the rate of production of cc pair integrated over the mass range 

9 Jlmr dMcc 

where is the charm quark mass and IniD is the DD threshold. 

The cross section for the low virtuality electroduction within CEM is 

^ep^e+Jl^^X ^ j^^^^ ^^2^1 dydxf^,Mhlp{^) (2.2) 

where fyj^iy) is the distribution function of the photon in the electron which, in the Weizsaker 
William approximation 1 13], is given by 

fy/eiy^E) = -{ 



ln{^-l\+\ 



71 y \ m 2 J 2 

{2-y)\ (2-2y 



}• (2-3) 



2y \2 — y 

To calculate SSA in scattering of electrons off a polarized proton target, we assume general- 
ization of this CEM expression for low virtuality electroproduction of //i/^by taking into account 
the transverse momentum dependence of the Weizsacker-Williams (WW) function and the gluon 
distribution function: 

^e+p^^e+J/w+X ^ t""^ dMlsdXydXg [d\i_yd\^g]fgip^{Xg,k^g) 
J Ami 

X/y/^(Xy,k_Ly) ■ (2.4) 

We assume dependence of pdf 's and WW function to be factorized in Gaussian form 

f{x,k^)=f{x)-^e-'^J^'^i) (2.5) 

with {k\) = 0.25 GeV^. 

The expression for the numerator of the asymmetry is 

d^r - d^r = 2L1 ^"^^ J \dxydxgd^k^yd'k^g]A'fgip,{xg,k^g) 

xf^/,{xy, k^y) 5\pg + py-q) al'^'\M^) (2.6) 
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where q = pc + Pc, A /^^^t {xg , kj^g) is the gluon Sivers function and M is invariant mass of the cc 
pair. 



The partonic cross section is [ |14j ] 

ao^^--"(M2) = \el^[il + r- ^r^)lni±^- (1 + 7)xA^] (2.1) 
where Y = 4ml/M^. 

Sivers asymmetry integrated over the azimuthal angle of J/y with a weight factor sin((^^j, — 
^s) is defined as 

/^^^r [/i"" [^M^]I[d^^±g]^'^fg/p^ {xgMg)fr/e{xr, qx - k^g)ao\sin{^g^ - ^s) 

An = (2.8) 

2 / d^qr U^^f [dM^] J [d^k^g]fg/p [xg , k^g)fy/e {xy, qr - k^g) do] 

where 0^^, and 05 are azimuthal angles of J/y and proton spin respectively and Xg^y = e^-^. 

3. MODELS FOR SIVERS FUNCTION 

We have used the following parameterization for the gluon Sivers function [^] 

Nf V./, A^^^\,,n. \ f . cos^k^. (3.1) 



A /g/pt {x,k±) = 2yKg (x) h{kx_ ) fg/p [x] 



There is no information available about the gluon Sivers function from experimental data. The 
valance and sea quark Sivers distribution functions used are the ones extracted from the HERMES 



and COMPASS experimental data in SIDIS processes [|15[]. 



The X dependent normalization for u and d quarks is given by, 

^,W^«,^'(l-.)'.^^i±M^ (.2, 

af fbf i 

where af,bf,Nf and Mi are best fit parameters obtained by fitting SIDIS, HERMES and COM- 
PASS data [§. 

For jyg{x), we have used two choices [ p^ ] 

(a) ^g{x) = {.jV,{^) + ^d{x))/2. 

(b) ^g{x)=^d{x). 

For h{k±^), we have used following two choices proposed by Anselmino etal ^: 

• Model I 

/.(^^) = ^/2;^.-^-^/^^ (3.3) 

• Modem 



where Mq = J {k"^) and Mi are best fit parameters. 
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Figure 1: (Color online) The single spin asymmetry A^'"^'''''^ "^^^ for the e + p^^e + J/\(/+Xasa function 
of y (left panel) and qj (right panel). The plots are for model I with parameterization (a) compared for 
JLab (V^ = 4.7 GeV) [soUd red line], HERMES = 7.2 GeV) [dashed green line], COMPASS = 
17.33 GeV) [dotted blue Une], eRHIC-1 (V? = 31.6 GeV) [long dashed pink Une] and eRHIC-2 (-^ = 
158.1 GeV) [dot-dashed black line]. 



4. NUMERICAL ESTIMATES 



We have used the following best fit parameters from the recent HERMES and COMPASS 



data [17 1 



A^„ = 0.40, fl„ = 0.35, bu = 2.6, 
Nd = -0.97, aa = 0.44, bd = 0.90 , 
M} = QA9GeV'^. 



(4.1) 



In figure 1 we have shown the comparison of y and qj distribution of estimated SSA at JLab, 
HERMES, COMRASS and eRHIC for model I and parameterization (a) of the gluon Sivers func- 
tion. The estimates are obtained using GRV98LO for gluon distribution function and Weizsaker- 
Williams function for photon distribution. The results for model II and parametrization (b) are 



given in reference 018|]. The hard scale involved in the calculation for all experiments is between 
Am^ and 4m|, as we are using color evaporation model. Hence the scale evolution of TMD's is not 
expected to affect much our estimates for the experiments at higher energies. 

According to our estimates sizable asymmetry is expected at various experiments covering 
different kinematical regions. Hence it is worthwhile to look at SSA's in charmonium production 
both from the point of view of comparing different models of charmonium production as well as 
comparing the different models of gluon Sivers function. 
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